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A novel technique for measuring the refractive index profile of optical 
fibers is demonstrated, which offers substantial advantages over alternative 
methods. The method consists of illuminating a small area of the fiber core 
and measuring the total transmitted power. The transmission of leaky 
modes is accounted for in the manner reported previously by other authors- 
The index profiles of germanium-doped fibers obtained by this technique are 
compared to interferometric measurements. The resolution is shown to be 
limited by wave optics effects to about \ {4n42A)~ l , where A = An/n. 
The distortion of the index profile as the wavelength varies and wave- 
optics effects are investigated. 

I. INTRODUCTION 

The accurate measurement of index profiles at various wavelengths 
may help design multimode fibers whose transmission capacity would 
go well beyond what has been presently achieved. Indeed, numerical 
calculations and theoretical analyses 1 - 2 show that there exist index pro- 
files (usually not power-law profiles), which, for quasi-monochromatic 
sources, provide transmission capacities of about 1.6/ A 2 Mb/s X km, 
where A = An/n. Measured transmission capacities are about 10 times 
smaller. To determine the optimum profiles, it is indispensable to know 
the variation of dn/d\ (where X is the operating wavelength) as a 
function of n for the class of materials considered with an accuracy 
of about 1 percent. The required variation of dn/d\ as a function of 
n can be obtained, in principle, from measurements on bulk samples 
(e.g., Ref. 3). We question, however, whether measurements on bulk 
samples are applicable to the fiber material with sufficient accuracy. 
For that reason and also because the fabrication and measurement of 
bulk samples is time-consuming, the direct measurement of index 
profiles at various wavelengths is highly desirable. Once the optimum 
profile applicable to the class of materials considered has been deter- 
mined, we measure the departures of the profile n(r) of the fabricated 
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fiber from optimum. Very small deviations may degrade considerably 
the transmission capacity. 

An interesting experimental technique for measuring circularly 
symmetric index profiles has been proposed by Gloge and Marcatili. 4 
The index profile is obtained by measuring with a pin hole the radial 
distribution of intensity in fibers excited by Lambertian (e.g., thermal) 
sources. In a series of careful measurements, Sladen, Payne, and 
Adams 8 have shown that good agreement can be obtained between 
the intensity in the fiber core and the index profile obtained by inter- 
ferometry provided the non-zero transmission of the leaky modes is 
accounted for. If this correction is made, the fiber samples need not be 
larger than about 1 meter, and may be as small as 1 cm. 

The technique that we describe in the present paper, which we 
call the transmission technique, is related to the near-field technique 
discussed above, but it differs from it in many significant ways. Arnaud 6 
has shown that, if we illuminate a small area of the fiber core (perhaps 
of the order of (\ 2 at x, y, the total transmitted power is, for sufficiently 
long fibers, proportional to n 2 (x, y) — n\, where n(x, y) denotes the 
refractive index at x, y, and n c , the cladding index. The proof is 
straightforward : The rays radiated from the illuminated area have an 
almost uniform distribution in the plane k z , k v , where k x , k y denote the 
transverse components of the wave vector k. Because of the relation 
kl + k 2 y + k\ = k 2 (x, y) = (o)/c) 2 n 2 (x, y), which holds between the 
rectangular components of k, and because only rays whose k z is larger 
than k s are transmitted without loss, the power transmitted through 
long fibers is proportional to 

k\ + k 2 = k 2 (x, y) - k\ = index profile (1) 

(see Fig. 1). The rays in the dotted area in Fig. lb leak away if the 
fiber is sufficiently long. Otherwise their contribution to the total 
transmitted power needs to be subtracted in the manner reported in 
Ref. 5. If the spot size is less than X , we may use as a, source either a 
(coherent) laser or a (spatially incoherent) led. If, however, the spot 
size is significantly larger than X , it is essential to use near-Lambertian 
sources such as leds. Indeed, coherent beams of large cross section 
would excite predominantly paraxial rays. This would require intro- 
ducing additional correction factors. 

In the present paper, we discuss the principles and limitations of 
the method, and we present experimental results. The transmission 
method gives results that are, in principle, identical to the near-field 
measurements described in Ref. 5. The main advantage of the trans- 
mission method, compared with alternative methods, including the 
near-field technique, is that it is extremely easy to implement. The 
results are highly reproducible to better than one part in 1000. 
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lb) 

Fig. 1 — (a) In the transmission method, the microscope objective illuminates a 
small area of the fiber end of the order of \ 2 . (b) The intensity is assumed uniform in 
the k x , k u space. For long fibers, only the rays in the shaded area are transmitted. 
For short fibers, the rays in the dotted area may also be transmitted (leaky rays). 
T denotes the power transmission. 

II. EXPERIMENTAL CONDITIONS 

To implement the proposed technique, all we need is an ordinary 
microscope, a high-radiance led (or a laser), and a detector. The 
numerical aperture (na) of the microscope objective should be at 
least twice as large as that of the fiber. One end of the fiber is properly 
broken or polished and centered approximately under the microscope 
objective at focus. When the microscope eyepiece is replaced by a led, 
a small spot of infrared radiation illuminates the fiber end. As we have 
indicated in the introduction, the power detected at the other end of 
the fiber is proportional to n 2 (x, y) - n 2 c , where n(x, y) denotes the 
index at the point x, y of the fiber where the light is focused, and n e 
the cladding index. To obtain the index profile, we may scan either the 
fiber, with a total motion of about 100 pm, or the source, with a total 
motion of about 3 mm. The arrangement shown in Fig. 2 incorporates 
a beam splitter (1) to allow the fiber to be observed during scanning. 
(A second beam splitter, which combines the light from two leds, is 
shown in Fig. 2. It is used only for dispersion measurements.) Some 
infrared leds radiate red light with sufficient intensity for direct visual 
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Fig. 2 — Experimental setup of the transmission technique. The fiber is scanned 
mechanically and its motion is recorded with a gauge. Two lbds are used for disper- 
sion measurement. 

observation. For leds at longer wavelengths, an image converter is 
necessary. Note that the focal length of microscope objectives may be 
slightly different for the red light and for the infrared light. In the 
arrangement in Fig. 2, it is convenient to have the distance between 
the led source and the microscope objective equal to the distance 
between the focal plane of the eyepiece and the microscope objective. 
This avoids the need for refocusing when the objective is changed 
from low to high magnification. To obtain good resolution, it is de- 
sirable that the led act as a point source; that is, that the apparent 
size of the led, demagnified by the microscope objective, be smaller 
than the diffraction-limited spot Fz5\ /na, defined by the numerical 
aperture of the microscope objective. An apparent emissive diameter 
of 25 nm (before demagnification) is adequate. The experimental set- 
up is shown in Fig. 3. The led source (not visible) is supported by the 
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xyz microscope stage at the left. The angular orientation of the fiber 
under the microscope objective can be varied. Figure 3 shows the 
differential micrometer that drives the fiber and the gauge that mea- 
sures its displacement with respect to the microscope objective. 

The advantages of the proposed technique compared to the more 
conventional near-field technique are many : 

(i) In the near-field method, the source is required to be Lam- 
bertian and uniform over the full cross section of the fiber core. 
As pointed out in Ref. 5, this condition is in fact difficult to 




Fig. 3 — Photograph of the experimental setup. The led (not seen) is supported by 
the xyz stage on the left. 
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achieve with leds. The authors in Ref . 5 thus found it necessary 
to use thermal sources instead of leds. Thermal sources (e.g., 
tungsten wires) provide poor signal-to-noise ratios when the 
spectral width is restricted by narrowband interference filters. 
In the transmission method, we require the na of the micro- 
scope objective to be significantly larger than that of the fiber 
(at least for coherent sources), but the requirement concerning 
the spatial uniformity of the source is relaxed. In some sense, 
the requirement of spatial uniformity is transferred from the 
source, where it is difficult to achieve, to the detector, where 
the condition is easily met. 
(it) The optics are much simplified. Only a single microscope 
objective is needed instead of (typically) three. Thus, the 
signal-to-noise ratio is improved. 

(Hi) Near-field measurements provide the shape of the index 
profile, but not the absolute value of An(r) = n(r) — n c , where 
n c denotes the cladding index. In the transmission method, we 
can calibrate An by measuring the intensity radiated axially 
by the microscope objective. This calibration technique will 
be discussed in more detail in the next section. 

(iv) The transmission method can be combined with the Fresnel- 
reflection technique (for a recent report of the Fresnel-reflection 
technique, see Ref. 7). To implement this modification, we 
replace the microscope eyepiece in Fig. 2 with a detector. 

An important drawback that applies to both the transmission and 
near-field methods is encountered when the fiber exhibits a low-index 
region near the cladding. In that case, some modes (besides the so- 
called weakly leaky modes) are leaking very slowly, and the inter- 
pretation of the measurements becomes ambiguous. The resolution 
offered by these methods may be marginal when the fiber profile 
exhibits very fast fluctuations. Note also that, for noncircularly 
symmetric profiles, the correction factors for leaky rays have not 
been worked out. If the deviation from perfect circular symmetry is 
small, however, the correction factor given in Ref. 5 may be used. 

III. INCIDENT BEAM PATTERN AND INDEX CALIBRATION 

To make precise measurements, the radiation from the microscope 
objective should approximately obey Lambert's law, at least for angles 
a to the axis that are less than V2A. To verify that this law is approxi- 
mately obeyed, we translate the detector in front of the microscope 
objective at some distance d J> \ from the focal point. Ideally, the 
variation of the detected power as a function of the distance r from 
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Fig. 4 — The curve shown is the desired radiation pattern from the microscope 
objective. It is sufficient if tliis law be obeyed from r/d = to the value correspond- 
ing to the A of the fiber (e.g., r/d < 0.3 if A = 0.02). The experimental points were 
measured using a Leitz UD 20, na = 0.57 microscope objective. 

axis should be* 

P(r) = cos 4 a = (1 + r 2 /d 2 )- 2 . (2) 

The desired variation of P with r in (2) is shown in Fig. 4. The maxi- 
mum value of r/d corresponding to a particular A is given by 

r/d = (na- 1 - I)"* (3a) 

na = n V2A. (3b) 

The values of r/d are shown in Fig. 4 for typical values of A and 
n = 1.45. 

Let us now consider the problem of calibrating An. This is done by 
measuring the intensity radiated axially by the microscope objective. 



* This result can be obtained by inverting eq. (5.246) of Ref. 6 and noticing that 
when the radiation is uniformly dense in the k x , k y plane (the transverse components 
of the wave vector), it remains uniformly dense after refraction at any plane interface 
perpendicular to the z axis. To show that, set in eq. (4.167) of Ref. 6 : dk x /dz = dk u /dz 
= (Descartes-Snell law) and find that/(A- I , k v , x, z) = g(k x , k v )h(x, z) is a solution 
of the Liouville equation for any differentiable functions g and h. Thus, if/ is inde- 
pendent of k z and k u at z = 0, as we have assumed (Lambert's law), it remains inde- 
pendent at kg, k„ after any number of refractions. 
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Let the power detected in front of the microscope objective be denoted 
P c and the power transmitted through the fiber for near-axial excitation 
be denoted P. If the detector radius is denoted p and its distance from 
the microscope objective focal point is d, the na = n V2A of the fiber 
is given by 

na = (j>/d)JFhP Cl (4) 

where 

v = [4n./(n. + l)»]» (5) 

accounts for the Fresnel reflection at both ends of the fiber. This ex- 
pression for 77 is not rigorous, but it is sufficiently accurate for our 
application. With sufficient accuracy, we can set n = 1.45. Then, 
rj = 0.93. Convenient values for d and p are d = 10 mm and p = 1 mm, 
respectively. We thus obtain from (4) and (5) 

na = 0.104VP/P C (6a) 

A = 0.00255 (P/P e ). (6b) 

It is, of course, necessary to have good breaks at both ends of the 
fiber. We have assumed that the fiber loss is negligible; this is the 
case for most fibers if the length is 1 m or less.* 

A more conventional technique for evaluating A consists of measur- 
ing the far-field pattern when the fiber is illuminated on or near axis. 
We have 

A = !(NA/n ) 2 , (7) 

where na denotes the sinus of the maximum radiation angle in air, 
defined typically at the — 3 dB point from maximum intensity. 

The core radius, a, of the fiber is best obtained by observing the 
fiber tip with a microscope. From the values of A, a, and X , the 
F-number is calculated according to 

V = (2to/Xo)na. (8) 

IV. CORRECTION FOR LEAKY RAYS 

If the fiber is not very long, the leaky modes excited by the source 
are not completely attenuated. They can be accounted for in the 
manner described in Ref. 5. Specifically, the index profile n 2 {r) — n 2 e is 
obtained by dividing the detected power P(r) by the correction factor 



* With some high na objectives (na ~ 0.8), spurious peaks are observed in the 
far-field pattern due to diffraction effects, even with led sources. Thus, the calibra- 
tion of A should be made with lower na objectives, e.g., na ~ 0.5. These spurious 
peaks do not appear to affect the profile measurements, but they prevented us from 
making a precise comparison between the two techniques described here for mea- 
suring A. 

1496 THE BELL SYSTEM TECHNICAL JOURNAL, DECEMBER 1976 



1 o 




X = 0.05/ 


f^OA 


1.5 


- 






■S 1-4 






0.15 


*-* 








o 








IT 

o 

1- 
u 

2 1-3 






^0.2 


z 
g 

O 
111 

gj 1.2 
O 

u 

1.1 

1 


- 


1 1 1 1 1 


>k 0.25 

0.5 

I 



0.2 0.4 0.6 0.8 1.0 

NORMALIZED RADIUS r/a 

Fig. 5 — Correction factor for leaky rays. 6 The X parameter is defined in the main 
text as X = V~ l log, (L/a). where V denotes the F-number of the fiber, L, the fiber 
length, and a, the core radius. The measured power should be divided by C. These 
curves are applicable to near-square-law fibers. 

C(r) given in Fig. 5. 6 In this figure, the parameter X is defined as 

X = V-nog e (L/a), (9) 

where L denotes the fiber length, and V the fiber F-number defined 
in (8). 

Note that the correction is negligible near the fiber axis, but may 
be as large as 30 percent at r Pi 0.8a for typical fiber lengths. Strictly 
speaking, the correction factor depends on the profile of the fiber. 
The curves in Fig. 5 are applicable to near square-law profiles. How- 
ever, the correction factor turns out not to be very sensitive to the 
exact profile. Thus, for most high-transmission capacity fibers, the 
square-law-profile correction factor may be sufficiently accurate. If 
greater accuracy is required, we may use an iteration procedure. 
Because this procedure is rather involved, it will not be discussed here. 

V. REFRACTIVE-INDEX PROFILE MEASUREMENTS 

The measurement technique described in previous sections has been 
applied to graded-index fibers. Let us first make a few general com- 
ments concerning the experimental technique and the results. The 
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results obtained are highly reproducible to better than one part in 
1000 even after a few hours if the fiber tip is protected by a glass 
plate and an index-matching fluid. The axial index dip characteristic 
of germanium- (or phosphor-oxide-) doped fibers is very useful to define 
the fiber center and achieve optimum focusing. 

5.1 Comparison with near-field and interferometric measurements 

Near-field measurements 5 were made on a graded-index fiber. The 
setup is shown schematically in Fig. 6a. A Burrus-type high-radiance 
led with an apparent emissive diameter of 50 /xm is imaged with unity 
magnification on the end of the fiber under test with a pair of micro- 
scope objectives (20 X, na = 0.4). The other end of the fiber is 
imaged with a microscope objective (40 X, na = 0.65) on a scanned 
small-area detector. The magnification, measured with a reticle, is 
equal to 54. Focusing is adjusted to make the details of the index 
profile as sharp as possible. The variation of detected power as a 
function of the transverse displacement of the detector is shown in 
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Fig. 6 — Near-field measurement on a germanium-doped, graded-index fiber. 
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Fig. 7 — Interferometric profile of the fiber in Fig. 6, X = 0.55 fim. na = 0.21, 
A = 0.0107. 8 

Fig. 6b for a fiber length L = 720 m. For such long lengths, the 
leaky-ray correction is small. Because of evaporation of germanium 
during the collapse of the preform, the material on axis is probably 
almost pure silica. The observed reduction of An is, at most, f . This 
discrepancy reflects the inherent limitation in resolution of near-field 
and transmission techniques. We attribute small oscillations on each 
sides of the central dip to wave-optics effects. They are qualitatively 
similar to the ones calculated for a dielectric slab in Appendix A. 
The slow modulation reflects the presence of diffused steps. The 
interferometric measurement 8 is shown in Fig. 7. 



OPTICAL FIBER INDEX PROFILE 1499 



L= 1.6 m 
X = 0.79 pm 



n 2 (r) - nj? (ARBITRARY UNITS) 



V\. 



-25 -20 



-5 

r IN jim 



10 



20 



25 



Fig. 8 — Profile of the fiber in Figs. 6 and 7 obtained with the transmission tech- 
nique. Uncorrected : solid line. Corrected for leaky rays : dashed line. Interferometric 
measurements from Fig. 7 : dots. The measured A at Xo = 0.79 /im is 0.0104. 

The index profile obtained with our transmission technique is 
shown in Fig. 8. The A of the fiber was measured from the far-field 
pattern, as described in Section III. We measured A = 0.0104, in good 
agreement with the value obtained by interferometry. From a core 
radius a = 24 jum, we calculate from (8) a V-number: V = 40. The 
fiber length is 1.6 m. Thus, the X-parameter in (9) is X = 0.28. The 
corrected profile, obtained by dividing the measured power by the 
correction factor C in Fig. 5, is shown as a dashed line in Fig. 8. The 
result of interferometric measurements is shown by dots for com- 
parison. Aside from the depth of the central dip, a significant difference 
between the dashed curve (corrected transmission profile) and the 
dots (interferometric measurement). Such a discrepancy may be in 
part attributed to the lack of perfect circular symmetry of the fiber. 

5.2 Measurement ot the lack of circular symmetry 

The (uncorrected) transmission profile of germanium-doped fiber 
was measured in two perpendicular azimuthal directions, labeled 0° 
and 90°, respectively. These two profiles are shown in Fig. 9 as solid 
lines and dashed lines, respectively. The measured A and na param- 
eters are given in the figure caption. We conjecture that, for the 
small deviations from circular symmetry exhibited by the fiber in- 
vestigated, the correction factor in Fig. 5 is applicable. 
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Fig. 9 — Uncorrected profile of a germanium-doped, graded-index fiber in two azi- 
muthal planes (0° and 90°). The measured na is 0.202, A = 0.00974, L = 1.6 m. 

5.3 Measurement of a double-hump profile fiber 

An attempt was made by MacChesney 9 to suppress the central 
index dip of germanium-doped fibers by increasing the amount of 
germanium halides during the final stage of fabrication of the preform. 
The profile of a fiber of that type was measured by interferometry 8 
and by our transmission technique. Both techniques clearly show that 
there is a large peak of index near the center of the fiber. However, 
the transmission technique shows that the dip in the center did not 
disappear (hence, the name "double-hump" given to the profile of 
that fiber). This central dip is not seen on the interferogram. The 
combination of a peak and a dip is unlikely to improve the fiber 
transmission. A much better compensation of the central dip will be 
reported later in this paper. 

5.4 Measurement of nonclrcularly symmetric profiles 

The theoretical result in (1) shows that the transmission technique 
is applicable, in principle, to noncircularly symmetric profiles, as well 
as to circularly symmetric profiles. A preform that accidentally 
collapsed flat 10 has been pulled at our request into a fiber and mea- 
sured. The uncorrected profiles are shown in Fig. 10. Because the 
correction for the leaky rays has not been made, the curves in Fig. 10 
are only indicative of the index profile. 
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Fig. 10 — Profile of a fiber with near-elliptical core (uncorrected). 

5.5 Profile distortion 

One of the most interesting and intriguing questions is whether 
index profiles get significantly distorted as the wavelength varies 
(independently of possible changes of scale). Fleming's measurements 
on bulk samples of germanium-doped silica 3 clearly indicate that pro- 
files get distorted significantly as the wavelength varies. This effect, 
however, has not been observed before on fibers. We report here 
preliminary measurements of profile distortion. The necessary for- 
mulas are relegated to Appendix B. 

Figure 1 1 shows the profiles measured with the transmission method 
at two wavelengths : X„ = 0.79 /xm and X = 0.9 /mi, on a germanium- 
doped, graded-index, large na fiber. Note first that the resolution 
(indicated by the depth of the central dip) is slightly poorer at the 
longer wavelength. When the scanning is made slightly off-center to 
avoid the central dip and the two profiles are normalized to unity on 
axis, the differences between the two profiles are found to be extremely 
small yet significant. To exhibit this difference with good accuracy, 
we have combined the fight from the two ledb with beam splitter 
in Fig. 2 number 2. Square pulses are applied to the ledb. The 
positive parts of the pulses drive one led and the negative parts drive 
the other. The levels are adjusted to have equal detected powers on 
the fiber axis, and therefore, zero signal on the lock-in amplifier. The 
difference between the two normalized profiles is plotted in Fig. 12 
(curve b). More precisely, we have plotted in Fig. 12 the "profile 
distortion" d = \ dr]/d\ , where n = N/2A and N m 1 - n 2 /n\, as a 
function of r/a. The accuracy of this curve is difficult to ascertain at 
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Fig. 11— Profile of a germanium-doped, graded-index fiber at two wavelengths (a) : 
X = 0.79 Mm, (b) X = 0.9 /mi, (c) X„ = 0.79, scanning slightly off-center (obtained 
with the transmission technique and an objective na = 0.85). 

the moment. Part of the observed change of the shape of the index 
profile may be attributed to the dependence of the leaky-ray correction 
factor on wavelength. 

The variation of d as a function of r/a can also be calculated from 
Fleming's measurements on bulk samples, as explained in detail in 
Appendix B. The result of this calculation is shown in Fig. 12 (curve a). 
There is no close agreement between curve a and curve c. However, 
they are comparable in magnitude. Thus, measurement of very slight 
changes of profile with wavelength, such as are shown in Fig. 12, are of 
great practical importance. 

VI. CONCLUSION AND PROPOSALS FOR FUTURE WORK 

We have proposed and demonstrated a novel technique for mea- 
surement of the index profile of multimode fibers. Conceptually, this 
technique is related to the near-field technique previously demon- 
strated by Sladen et al. 5 (near-field technique). From a practical 
point of view, however, our technique is quite different, since it does 
not require Lambertian sources. In particular, lasers can be used. We 
have found that the measured profiles are highly reproducible, to 
better than one part in 1000 over periods of hours. Index profile 
measurements can be obtained in a few minutes including fiber-end 
preparation. The agreement between our technique and interferomet- 
ric measurements leaves something to be desired. The discrepancy, 
however, may be attributed to the lack of perfect circular symmetry 
of the fibers investigated. Theoretical considerations show that the res- 
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Fig. 12 — Profile distortion. Curve (a) is calculated from Fleming's measurements 
on bulk samples as reported in Ref. 3 (Fig. 1, curve labeled X = 0.9 jum). Curve (b) 
is the difference between the normalized profiles at Xo = 0.79 and 0.9 /*m. Curve (c) 
is the same as curve (b) corrected for the leaky rays. 



olution is about X/4V2A. For a typical value A tt 0.015, this resolu- 
tion is about the free-space wavelength X„ ftJ 1 pm. This appears to be 
sufficient for most practical purposes. 

Comparison of depths of central-index dips suggests that the trans- 
mission technique (and the near-field technique as well) provides 
better resolution than interferometric techniques. We have pre- 
sented preliminary evidence for the distortion of the index profile as 
the wavelength varies (profile distortion), an effect that was inferred 
previously only from measurements on bulk samples. Theories that 
neglect profile distortion may be in considerable error. 

We shall now make a few suggestions for improvement of the 
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measurement technique. Oil-immersed microscope objectives would 
be useful to prevent interference effects between the objective and the 
fiber tip when monochromatic laser sources are used. The processing 
of the experimental data can be considerably improved if the lock-in 
amplifier and the gauge have digital readouts. The correction factor 
C, given approximately in Fig. 5, should be recalculated and supplied 
in digital form. A better approach would consist of calculating C from 
the apparent measured profile and iterating. These iterations, pre- 
sumably, can be effected with modest computing time. Noncircularly 
symmetric profiles can be corrected, in principle, but the correction 
problem has not been solved yet. Finally, we attempt to deconvolve 
the wave-optics effects (discussed in Appendix A) that are most 
conspicuous in regions where the index varies rapidly. The possibility 
of performing this deconvolution is intriguing, but the analytical 
problem remains, to our knowledge, unresolved. The case of fibers 
with an index barrier between the core and the cladding requires 
further analysis. 

Among all the index-profile measurement techniques that have been 
proposed so far, the transmission technique that we have described 
here appears to be the easiest to implement and the most reliable. 
Improvement in data processing should make the results quite accurate 
in most cases. 

APPENDIX A 
Wave Optics Effects 

The result (given in the main text) that the intensity distribution 
in the cross section of long multimode fibers is proportional to 
k 2 (x, y) — k 2 s , where k(x, y) m (w/c)n(x, y) denotes the core wave 
number and k, = (co/c)n c the cladding wave number, is based on ray 
optics (wkb approximation) and on the assumption that rays whose 
axial wave number (k z ) is less than k, are radiating away and do not 
contribute to the total intensity. Because the number of trapped modes 
carried by real fibers is finite, the intensity distribution does not follow 
the fine details of the index profile. This is because the optical field 
cannot vary in transverse directions faster than (sin k x max :c) where 
jc im „ = (a;/c)n V2A denotes the maximum value of the transverse 
wave number, n the index on axis, and A £d An/n. According to the 
above formula, the smallest distance between nodes and peaks of the 
irradiance in the fiber core is 

Ax«X„/(4n V2AJ. (10) 

Equation (10) provides an estimate of the resolution afforded by the 
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method. For example, if X = 1 /mi, n = 1.45, and A = 0.01, the reso- 
lution is, according to (10), Ax tt 1.25 /mi. 

To get a more precise estimate of the error at discontinuities of the 
index profile, let us consider, as a model, an oversized dielectric slab. 
When the origin of the x axis is taken at one of the slab boundaries, 
the normalized field of //-modes of order m is (see, for example, Ref. 
6 with a slight change of notation) 



E m {x) = i 
where 



- ^ sin (Y m X - arcsin Y m ), X £ 0, 
^exp(-VT^nX), Z^O, 



(11) 



Y m = (m+l)(7r/2)F- 1 (12) 
and X, V are normalized distances and frequencies 

X = /cV2Ax (13a) 

V = fcV2Ad. (13b) 
For trapped modes, we require 

Y m < 1. (14) 

When the fiber is excited by a Lambertian source, the intensity I(x) 
in the slab cross section is obtained by adding the intensities, E^(x), 
of all the trapped modes. Because the slab is highly oversized, we can 
replace the summation over m by an integral. We obtain 



■ir/2 

I(x) = 



\2 r sin 2 (sin dX - 6) cos Odd X ^ 

J ° (15) 

/•t/2 V ' 

2 / sin 2 6 exp (-2 cos 6X) cos ddd, X ^ 0. 
- Jo 



For X — ► — °o and X — > + 00 , we have as expected, 1 = 1 and 1 = 0. 
At the discontinuity (X = 0), an elementary integration of (15) gives 
/ = f . The intensity profile defined in (15) is shown in Fig. 13. We 
notice an overshoot of the irradiance equal to 14 percent. This over- 
shoot is somewhat similar to the Gibbs effect encountered with Fourier 
series. If we keep the slab thickness a constant but increase the optical 
frequency, the region where the irradiance departs significantly from 
the ray-optics value becomes narrower and narrower. The amplitude 
of the overshoot, however, remains the same. The curve in Fig. 13 
provides understanding of the limitation in resolution of the method 
discussed in the main text. This limitation is caused by wave-optics 
effects that have been ignored in eq. (1). 
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Fig. 13 — Wave-optics effect in highly multimoded slab (a). The curve in (b) 
exhibits the departure from ray optics. Note the overshoot of 14 percent, which is 
independent of wavelength. 

APPENDIX B 
Profile Distortion 

In this Appendix we derive formulas relating to the change of shape 
of the index profile of a fiber as the wavelength varies (independently 
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of a possible change of scale). These formulas relate the distortions 
measured directly on the fiber to measurements on bulk samples. 
Let us assume that we have normalized the profiles to unity at r = 0. 
The ordinate of the curve is 1 — 17, where 77 = N/2A, N(r) = 1 
— n 2 (r)/n 2 and 2 A = 1 — n 2 c /n 2 . From the normalized profiles at two 
closely spaced wavelengths, we evaluate the profile distortion d(r) 

d(r) = \ o dri/d\ . (16a) 

This may be written as 

d = 0" l l\o(dP/d\„) + P(X/C)(dC/dX), (16b) 

where P(r, \ ) is the normalized (but uncorrected) Siber transmission, 
C(r, X) the leaky-ray correction factor, and X is the parameter 
defined in (9). 

C and the differential correction term (X/C)(dC/dX) can be 
obtained from Fig. 5. 

Let us now consider the curve S(n) where S =—\ ndn/d\ that 
can be obtained from measurements on bulk samples. From the value 
of A and the cladding index n c , we calculate the index on axis n . It 
is not difficult to show [e.g., from eq. (14) of Ref. 11 ] that the dis- 
tortion parameter d defined in (16) is related to the dispersion param- 
eter S defined above by 

d( v ) = lD o /A(n 2 + S )Ts(v) ~ *(D1 (17) 

where we have defined 

D m 1 - (\ /n )(dn o /d\ o ), (18) 

and 

8 = S(n) - S(n ) (19) 

is considered a function of 77. In particular, s(l) = S(n c ) — S(n ). It 
is clear, from its definition and from (17), that 

d(0) = d(l) = 0. (20) 
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